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ABSTRACT 


The  equation  defining  the  energy  radiated  per  unit  solid  angle  due  to  Cerenkov 
radiation  is  analyzed  in  detail,  including  the  effects  of  all  equation  variables  for  a 
hypothetical  electron  accelerator  experiment.  Specifically,  various  finite  frequency 
intervals  are  used  in  an  effort  to  determine  the  optimum  means  of  determining  the 
details  of  a  charge  bunch  in  a  high  energy  electron  accelerator.  In  particular,  it  is 
shown  how  narrowband  measurements  as  a  function  of  angle  may  yield  information  on 
both  the  beam  path  length  and  the  bunch  charge  parameters.  As  an  aid  to  the 
analysis,  an  interactive  Fortran  program  is  presented  that  allows  for  any  specific 
experimental  parameters,  with  options  for  various  output  types  as  desired. 


THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may 
not  have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made, 
within  the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and 
logic  errors,  they  cannot  be  considered  validated.  Any  application  of  these  programs 
without  additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


A.  HISTORY 

The  existence  of  the  phenomenon  now  known  as  Cerenkov  radiation  was 
observed  as  early  as  1910.  most  notably  by  Madam  Curie.  However,  during  the 
ensuing  years,  there  was  other  work  being  completed  by  those  familiar  with  Madam 
Curie's  observation  that  masked  a  detailed  study  of  the  phenomenon.  Jelley  [Ref.  1] 
explains  their  reasons  in  greater  detail.  In  1926.  Mallet  made  the  first  deliberate 
attempt  to  study  and  explain  “the  very  faint  emission  of  a  bluish-white  light  from 
transparent  substances".  Experimentation  concerning  this  phenomenon  continued 
through  the  I920's  and  1930  s  by  Cerenkov  and  Mallet,  but  a  viable  theory  explaining 
the  process  would  not  be  proposed  until  1937  by  Frank  and  Tamm.  Their  theory,  for 
which  they  were  consequently  awarded  the  1958  Nobel  Prize  in  Physics,  was  found  to 
be  in  excellent  agreement  with  the  experimental  results  obtained  by  Cerenkov.  The 
advent  of  more  sensitive  light  detectors  and  other  experimental  equipment  accounted 
for  more  in  depth  studies  of  Cerenkov  radiation  in  the  1940's  and  1950  s.  [Ref.  1:  pp. 
1-8] 

Since  then,  numerous  experiments  have  been  conducted  and  many  papers  and 
academic  theses  published  that  have  helped  to  more  clearly  understand  and  explain  the 
Cerenkov  radiation  phenomenon.  Some  of  the  more  recent  investigations  include; 
periodic  electron  bunches  of  finite  emission  lengths,  time  developement  of  bunch 
charges  and  Cerenkov  radiation  in  the  x-ray  region,  (see  Bibliography). 

B.  BACKGROUND 

Work  on  this  thesis  was  motivated  by  two  factors.  First,  using  previously  derived 
equations  and  expressions  describing  the  energy  radiated  per  unit  solid  angle  due  to 
Cerenkov  radiation,  [Ref.  2]  a  fortran  program  was  written  to  run  on  the  IBM  3033 
mainframe.  The  program,  (see  Appendix),  was  written  interactively  to  allow  all 
parameters  to  be  changed  in  order  to  tailor  the  output  to  fit  any  specific  experimental 
setup. 

Second,  using  the  aforementioned  Fortran  program,  a  series  of  outputs  were 
generated  to  be  used  in  comparing  the  theoretical  results  of  previously  derived 
expressions  with  those  of  actual  experimentation.  Since  there  are  literally  an  infinite 


II.  TIIEORV 


A.  CERENKOV  EEFECT 

Cerenkov  radiation  can  most  easily  he  explained  by  describing  the  interaction  of 
a  single  electron  within  a  transparent  medium.  Consider  a  mm-relativistic  electron. 
'■Chile  traversing  the  medium,  the  electron  will  tend  to  locally  polarize  the  adjacent 
atoms  m  the  medium,  instantaneously  creating  a  temporary  dipole  within  tire  material. 

I  i icrc: cro  as  the  electron  moves  through  the  medium,  an  electromagnetic  pulse  is 
generated.  However,  since  there  is  complete  symmetry,  there  will  he  no  net  field 
generated.  See  I  .cure  2  1 


(A)  IB) 

Figure  2.1  Polarized  Atoms  in  a  Dielectric. 

Now  suppose  an  electron  is  travelling  at  a  speed  comparable  to  or  exceeding  the 
speed  of  light  lor  the  medium.  Although  symmetry  is  preserved  in  the  azimuthal  plane, 
along  the  axis  of  motion  a  resultant  dipole  field  is  generated,  these  fields  are  set  up  at 
each  element  along  the  electron  s  track,  radiating  a  brief  electromagnetic  pulse. 
Provided  the  electron  is  moving  at  a  speed  greater  than  the  speed  of  light  within  the 
medium,  the  wavelets  formed  from  all  elemental  positions  on  the  track  can  be  in  phase 
and  thus  produce  a  resultant  electromagnetic  field.  1  Ref.  I:  pp  ?-b| 

1  igure  2  2  depicts  the  relationship  ol  the  distance  iHZ)  travelled  by  an  electron  in 
time  At  versus  the  distance  iB  1)  covered  by  an  emitted  electromagnetic  pulse  during 
the  same  time.  1  he  distance  ( DZ )  travelled  by  the  electron  is  given  bv; 


BZ  =  Pc,, At. 


(cqn  2.1) 


53 
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where  j)  is  the  ratio  of  the  electron  speed  to  the  speed  of  light  in  a  vacuum  ,cQ.  The 
distance  (/?.!)  travelled  by  the  electromagnetic  wave  is  given  by: 


BA  =  cAt, 


(eqn  2.2) 


where  c  is  the  speed  of  light  for  the  medium.  For  a  medium  with  an  index  of  refraction 
n,  c  =  c0/n.  From  equations  2.1  and  2.2,  the  “cerenkov  relation"  is  obtained. 
[Ref.  I:  pp  6-S| 


cos0c  *  1/pn. 


VAt 


(eqn  2.3) 
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Tigi^re  2.2  Illustration  of  Cerenkov  Radiation. 

Due  to  the  finite  path  length  at  which  the  radiation  is  detected,  the  Cerenkov 
cone  angle  is  normally  shifted  away  from  0C.  As  the  frequency  at  which  the  emitted 
radiation  is  measured  or  the  distance  from  the  source  is  increased,  the  radiation  cone 
angle  approaches  0c.  As  this  occurs,  an  increasing  fraction  of  the  total  Cerenkov 
radiation  is  found  at  0c.  [Ref.  2J 


WJWt 


B.  MATHEMATICAL  INTERPRETATION 
1.  Energy  Equation 

As  previously  discussed,  for  a  finite  path  length,  diffracted  Cerenkov  radiation 
effects  will  be  seen  at  angles  other  than  the  Cerenkov  angle.  Previous  work  on  this 
phenomenon  (Ref.  2]  has  resulted  in  an  expression  for  the  emitted  Cerenkov  radiation. 
The  energy  radiated  per  unit  solid  angle  within  the  frequency  range  dv  by  a  single 
bunch  charge  q  travelling  a  finite  distance  L  is: 

E(v,X:)dv  =  Q R*dv,  (eqn  2.4) 

where  Q  is  a  constant  whose  magnitude  is  dependent  on  the  index  of  refraction  and  the 
total  charge  in  Coulombs. 


Q  =  \lcq2/8K2.  (eqn  2.5) 

The  radiation  function  R  is  given  by: 

R  =  27tqsin0l(u)F(£),  (eqn  2.6) 

where  0  is  the  angle  measured  from  the  line  of  travel  of  the  bunch  charge  to  the 
direction  of  propagation  of  the  emitted  radiation,  I(u)  is  the  diffraction  function,  and 
F(  A )  is  a  form  factor.  For  convenience,  the  length  of  travel  of  the  bunch  charge  is 
measured  in  units  of  the  wavelength  of  the  emitted  radiation  within  the  medium.  The 
dimensionless  beam  length  parameter  q  serves  this  purpose  and  is  defined  as  the  ratio 
of  the  length  L  from  the  source  to  the  wavelength  X.  of  the  emitted  radiation. 

q  =  L/X.  (eqn  2.7) 

The  diffraction  function  is  defined  as: 

l(u>  =  sin(u)/u,  (eqn  2.8) 


where  u  is  dependent  on  both  the  angle  0  and  the  beam  length  parameter  q. 


/1/^r 


Equations  2.5  through  2.9  can  then  be  used  to  redefine  equation  2.4  in  terms  of 
experimental  constants,  (i.e.  n,  L,  P  etc.)  and  the  single  remaining  variable  of 
frequency.  [Ref.  2] 


E(v,A:)dv  =  QG2F2(£)sin2(Av}dv. 


(eqn  2.10) 


For  a  given  experimental  setup,  A  and  G  are  functions  of  the  angle  only  and  are 


defined  as  follows: 


A  =  {L7t/c}*{(l/nP)-cos8},  and 


(eqn  2.1 1) 


G  =  {2sin0}/{(l/np)-cos0}. 


(eqn  2.12) 


The  radiation  patterns  arising  from  equation  2.4  can  be  thought  of  as  an 
oscillating  sin2x  function  modulated  by  an  envelope.  If  F{k)  is  neglected,  then  G2  acts 
as  the  modulating  envelope.  The  actual  form  of  the  envelope  will  vary  with  nP,  but 
will  be  constant  for  a  given  experimental  setup.  Figure  2.3  show's  a  plot  of  G2  as  an 
illustrative  example  for  values  of  nP  two  percent  above  and  below'  threshold, 
(nP=  l.0±0.02). 

Finally,  the  dimensionless  form  factor  F(/c)  is  related  to  the  Fourier  transform 
of  the  bunch  charge  distribution.  The  Fourier  components  of  any  bunch  charge  are 


defined  as: 


p(/c)  *  Jp^e^dV, 


(eqn  2.13) 


where  p (r)  defines  the  charge  distribution  of  the  bunch  charge.  Once  the  Fourier 
components  are  determined  from  equation  2.13,  the  form  factor  is  found  by  the 


relation: 


p(/c)  =  qF(/c). 


(eqn  2.14) 


Here  we  consider  only  a  line  charge  with  distribution  p(z)  which  is  travelling  in  the  z 
direction.  Thus,  r-*z  in  equation  2.13  and  &-*kz  in  equation  2.14.  In  this  work,  only 


fWAV.V.V 


the  trapezoidal  bunch  charge  distribution  is  used  although  many  others  are  available 
and  easily  used  within  the  Fortran  program.  [Ref.  3] 

2  Trapezoidal  Form  Factor 

The  first  step  in  obtaining  any  type  of  form  factor  is  to  determine  the  assumed 
geometry  of  the  emitted  bunch  charge.  Figure  2.4  depicts  the  geometry  for  the 
trapezoidal  bunch  charge.  The  functions  defining  the  positive  half  of  the  trapezoid  are 


given  by: 

f,(x)  =  A,  and  (eqn  2.15) 


f2(x)  =  (A/(d-b))  •{d-x> ,  (eqn  2.16) 

where  fj  defines  the  trapezoid  from  0  to  b  and  f2  defines  the  trapezoid  from  b  to  d. 
The  amplitude  of  the  bunch  charge  is  denoted  as  A  with  units  of  coulombs  per  meter. 
The  values  b  and  d  have  units  of  meters  and  arc  easily  obtained  from  the  pulse  length. 


Figure  2.4  Trapezoidal  Bunch  Charge. 

The  fourier  components  of  this  trapezoid  arc  defined  by  equation  2.13  where 
ptk^)  is  the  sum  of  the  two  functions  f((x)  and("2(x),  (ie.  the  charge  distribution  of  a 
single  bunch  charge),  and  k2  is  the  component  of  the  wave  vector  k.  The  details  of 
solving  equation  2.13,  though  tedious,  arc  rather  straightforward.  [Ref.  3) 


Equation  2.17  is  the  solution  to  equation  2.13  and  provides  the  Fourier  components 
in  the  case  of  the  trapezoid. 


p(kz)  =  (  (2A)(d-b)/k*z  )•{  cos(kzbi  -  cos(kzd)  }.  (eqn2.l'T) 

The  value  of  A  is  found  by  equating  q,  the  total  bunch  charge,  to  the  integral 
of  the  charge  distribution. 

q  =  _  xjX  P(z>dz  =  A(b  +  d).  (eqn  2. IS) 

By  substituting  the  results  of  equations  2.17  and  2.18  into  equation  2.14.  the  form 
factor  for  the  trapezoidal  charge  distribution  is: 

F(kz)  =  {  2q/(d2-b2)  } •{  1/k2) •{  cos(kzb)-cos(kzd)  }.  (eqn  2.19) 

If  the  medium  is  assumed  to  be  non-dispersive,  then  to  =  ckz  and  remembering 
that  kz=  AcosB,  the  form  factor  can  be  written  in  terms  of  the  angular  frequency 
2ttv  =  o).  After  substituting,  the  explicit  relation  for  the  trapezoidal  form  factor  is 
given  by  equation  2.20  as  follows: 

F(kz)  =  {  c2,2n2(d2-b2)  }•{  cos(Bv)-cos(Dv)  }•(  1/v2  },  (eqn  2.20) 
where  B  and  D  are  given  by  equations  2.21  and  2.22  respectively. 

B  =  2rcbcos0/c,  and  (eqn  2.21) 


D  *  2ttdcos0/c.  (eqn  2.22) 

It  is  clearly  seen  that  once  the  parameters  for  the  bunch  charge  are  defined, 
the  form  factor  is  highly  dependent  on  the  inverse  square  of  the  frequency  v.  The 
expression  for  the  energy  radiated  per  unit  solid  angle  within  the  frequency  range  dv 
written  in  terms  of  frequency  will  then  be  equation  2.4  with  the  form  factor  as  given  in 
equation  2.20.  [Ref.  4] 
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(eqn  2.23) 
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Edv  =  [  QG2sin2(Av)  ] 

•[  {c2/2n2(d2-b2)}  (cos(Bv)-cos(Dv)}  (1/v2}  ]2dv. 

The  radiated  energy  is  a  function  of  frequency  through  both  the  sur(Av)  term  and  the 
form  factor. 
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III.  RESULTS  AND  ANALYSIS 


The  Fortran  program  written  for  this  work  has  various  output  capabilities.  The 
first  type  of  output  contains  data  equating  the  explicit  value  of  E  (eqn.2.4),  in  joule- 
seconds.  to  a  specific  frequency  at  a  particular  angle  8.  Another  output  type  provides 
data  for  graphing  the  radiation  envelope  (eqn.2.3)  or  the  form  factor  (eqn.  2.20)  used  in 
this  thesis.  These  graphs  are  useful  in  understanding  the  complex  nature  of  the 
variables  involved  in  solving  equation  2.4. 

The  most  useful  output  type  is  that  which  solves,  (i.e.  integrates),  equation  2.4 
over  a  predetermined  frequency  range  for  all  desired  angles  from  the  beam  path.  A 


careful  study  of  these  graphical  outputs  will  provide  insight  into  the  particulars 
necessary  for  any  experimental  setup  so  that  the  most  accurate  and  informative  results 
can  be  obtained. 

As  previously  discussed,  the  output  obtained  during  this  work  are  based  on  a  set 
of  parameters  that  uniquely  define  a  specific  electron  accelerator  experimental  setup. 
The  following  input  parameters  (with  their  Fortran  variable  names),  are  required: 

1.  Accelerator  Beam  Energy  (EBEAM). 

2.  Total  Periodic  Electron  Bunch  Charge  (CUE). 

3.  Index  of  Refraction  (IND). 

4.  Distance  from  Source  to  Detector  (LENGTH). 

5.  Cerenkov  Radiation  Frequency  Interval  (MIDNU.ENU). 

6.  Incremental  Change  in  Frequency  Interval  (DNU). 

7.  Angle  of  Interest,  as  Measured  from  the  Charge  Path  (DTHETA). 

8.  Bunch  Charge  Form  Factor  (FORFAC). 

9.  Trapezoidal  Bunch  Parameters  (BEE, DEE). 

The  specific  values  for  each  variable  used  to  obtain  the  outputs  in  this  work  are 
as  follows: 

1.  Beam  Energy  -  50  ME V. 

2.  Bunch  Charge  -  .001  Coulombs. 

3.  Index  of  Refraction  -  1.000268  (air). 


Source  to  Detector  Distance  -  100  Meters. 


Frequency  Interval  -  As  noted. 
Frequency  Increment  -  0.5  MHz. 


7.  Solid  Angle  -  As  noted. 

8.  Form  Factor  -  Unity  and  Trapezoid. 

9.  Trapezoidal  Bunch  Parameters. 

a.  Top-50  nanoseconds. 

b.  Base-60  nanoseconds. 

For  the  given  index  of  refraction  and  beam  energy,  the  “Cerenkov  Relation"  (eqn.2.3) 
gives  the  critical  Cerenkov  angle  as  1.9°. 

Three  basic  assumptions  were  made  in  the  course  of  this  work.  First,  it  is 
assumed  that  the  shape  of  the  trapezoidal  bunch  charge  remains  unchanged  as  it 
travels  through  the  medium.  This  is  reasonable  since  the  distances  involved  are 
relatively  short  with  respect  to  the  speed  at  which  the  bunch  charge  is  travelling. 
Second,  the  permeability  of  free  space  p  is  used  vice  that  of  the  actual  medium  of  air. 

The  final  assumption  is  that  the  speed  of  the  electron  bunch  remains  unchanged 
throughout  it's  travel.  Although  there  is  some  Bremstrahlung  radiation  emitted,  it  has 
a  negligeable  effect  on  the  speed  of  the  bunch  charge  and  therefore  P  is  assumed  to 
remain  constant. 

Figure  3.1  depicts  the  energy  radiated  per  unit  frequency  at  angles  of  30.  45  and 
60  degrees  using  a  form  factor  of  unity.  As  would  be  expected  from  equation  2.23,  a 
sin2x  function  results  since  Q  (eqn.2.5),  A(eqn.2.ll)  and  G(eqn.2.12)  are  all  constant  at 
the  given  angle  0,  the  only  variable  being  the  frequency  v.  The  difference  in 
amplitudes  is  due  to  the  variation  of  G  with  0.  The  difference  in  the  periods  is  due  to 
the  function  A,  which,  as  the  argument  of  the  sin  term,  effects  the  periodicity. 

Once  a  form  factor  other  than  unity  is  used,  the  output  pattern  is  changed 
significantly.  Equation  2.23  shows  the  frequency  dependence  of  the  trapezoidal  form 
factor.  Figure  3.2  is  the  form  factor  obtained  at  0  degrees  over  the  frequency  range 
10-100  MHz.  The  damping  of  the  form  factor  waveform  is  due  to  the  inverse  square 
relationship  of  F (k)  with  frequency.  Figures  3.3  and  3.4  depict  the  square  of  the  form 
factor  as  it  is  used  in  equation  2.4.  Different  frequency  intervals  are  used  but  each 
waveform  is  computed  at  the  same  angle,  0  =  0°.  As  with  Figure  3.2,  Figures  3.3  and 
3.4  each  dampen  with  increasing  frequency.  The  increased  periodicity  of  Figure  3.4  is 
due  primarily  to  the  wider  frequency  range.  At  angles  other  than  0  =  0°  the  form 
factor  will  be  similar  in  shape  except  scaled  by  changes  in  B  and  D  due  to  the  cos0 
factor.  [Ref.  4] 


ENERGY  DISTRIBUTION  (F(K)=1.0) 


f  igure  3.1  Energy  Radiated  per  Unit  Frequency,  (0  =  30°, 45°  and  60°) 


Figure  3.5  is  a  plot  similar  to  that  of  Figure  3.1  except  that  the  unity  form  factor 
has  been  replaced  with  that  of  the  trapezoid.  As  expected,  the  use  of  a  specific  form 
factor  significantly  effects  the  output.  The  output  of  interest  wii!  be  the  integration  c f 
ail  outputs  similar  to  Figure  3.5  at  each  angle  of  interest.  Beiore  anaiv/me  the 
integrated  output  with  the  trapezoidal  form  factor,  it  useiul  to  look  at  1  mures  3  o 
through  3.S  which  depict  the  integration  of  equation  2.2?  while  maintaining  units  Arm 
factor.  As  noted  earlier,  an  energy  peak,  at  the  Cerenkov  angle  of  l.\T  is  expected. 
Although  the  peak  appears  to  be  centered  at  about  In  '  m  Figure  ?.<■.  as  the  frequency 
interval  increases  from  lo-luo  MHz.  to  10- io.onn  MHz.  the  same  curse  shifts  closer  to 
the  Cerenkov  angle.  Concurrently,  the  graphs  become  smoother  as  more  data  points 
are  used  at  each  interval  0.  The  aberration  found  in  Figure  3 A  may  be  an  artifact  o: 
the  Fortran  program.  The  extraordinarily  large  magnitude  of  Figure  3.5  is  due  to  the 
existence  of  a  pole  m  variable  G  ieqn.2.12)  of  equation  2.10  as  0  approaches  the 
critical  Cerenkov  angle.  However,  since  the  variable  A  is  zero  at  the  Cerenkov  Angle, 
equation  2.10  remains  finite.  [Ref.  2j 

Figure  3.9  shows  the  output  of  equation  2.23  integrated  over  a  relatively 
broadband  frequency  range  of  100  to  1.000  MHz.  Figure  3.10  is  the  identical  output 
plotted  over  a  narrower  range  of  8  to  enhance  the  output  detail.  In  each  case,  the 
detail  of  the  output  is  fuzzy  at  best,  making  it  difficult  to  conduct  any  meaningful 
analysis.  In  an  attempt  to  force  more  detail  in  the  output,  the  frequency  range  over 
which  equation  2.23  is  integrated  was  narrowed.  Figure  3.11  depicts  the  energy  spectra 
integrated  ever  20  MHz.  centered  around  100  MHz.  All  other  parameters  are  equal  to 
those  used  in  Figures  3.9  and  3.10.  Although  similar  in  general  form  .  Figure  3.11 
shows  more  detail.  Then,  changing  to  a  narrower  frequency  band  of  2  MHz. (again 
centered  around  100  MHz.),  Figure  3.12  depicts  a  clearly  defined  output  that  is  similar 
to  a  frequency  modulated  carrier  wave.  Figure  3.13  shows  the  same  output  lor  a 
narrower  range  of0  to  show  more  detail. 

The  Cerenkov  radiation  envelope  discussed  earlier!  Figure  2.3i  has  an  impact  on 
the  outputs  of  Figures  3.12  and  3.13,  but  once  nP  is  selected  the  envelope  will  effect 
only  the  amplitude  of  this  and  similar  graphs.  The  null  points  of  Figures  3  12  and  ?  13 
can  then  be  calculated  by  analyzing  the  components  of  the  radiation  function!  eqn. 2. m. 
The  expression  for  the  trapezoidal  bunch  charge  form  factor!  eqn.  2. 20)  can  be  expressed 
as  a  product  of  a  constant  and  two  sin  functions  using  the  following  trigonometric 
identity. 
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ENERGY  DISTRIBUTION/!  (F(K)=  TRAPEZOID) 


cosa  -  cosy  =  -2sin{  (a,2)  +  (y  2)  }*sin{  (a  2)-(y -2)  }. 


(eqn  3. 1 ) 


After  a  great  deal  of  algebraic  manipulation,  it  can  be  shown  that  the  form  factor  will 
be  of  the  form: 

F(A)  =  const*. sin(£)*.sin(<p),  (eqn  3.2) 

where  c  =  rtcos0(d+ b)v  c.  and  (p  =  rtcos0(d-b)v  c. 

For  the  trapezoidal  bunch  charge  parameters  used  in  this  work,  c  (d-b)  equates 
to  a  value  of  200  MHz.  and  c  (d-t-b)  to  a  value  of  18.2  MHz.  Then,  using  the  mean 
frequency  of  100  MHz.  found  in  Figure  3.13,  one  need  only  to  equate  each  of  the  sin 
function  arguments  to  integer  values  of  rr  to  mathematically  determine  the  values  of  0 
at  which  equation  3.2  will  equal  zero.  Table  1  lists  the  values  of  0  at  which  the  form 
factor  will  force  zeroes  in  the  energy  equation. 


TABLE  1 

VALUES  OF  0  AT  WHICH  THE  FORM  FACTOR  EQUALS  ZERO. 


sinc^O 

sin(p  =  0 

0j  =  79.5°. 

undefined. 

02  =  68.7°. 

undefined. 

03  =  56.9°. 

undefined. 

04  =  43.3°. 

undefined. 

05=  24.5°. 

undefined. 

At  the  mean  frequency  of  100  MHz.,  sin  <p  =  0  cannot  be  determined  since  it 
requires  solving  for  an  angle  0  whose  cos  is  greater  than  unity.  However,  by 
comparing  the  values  of©  in  column  1  of  Table  1  to  Figure  3.13,  one  sees  a  correlation 
of  the  null  points  to  what  can  be  considered  the  envelope  of  the  output.  It  is  clear 
then  that  the  form  factor  acts  as  the  modulator  of  the  energy  spectra  of  the  Cerenkov 
radiation. 
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figure  3.9  Energy  per  Unit  Solid  Angle,  lOO-IOOOMHz.,  (TrapczoiJal  Form  Factor) 


TABLE  3 

VALUES  OF  0  AT  WHICH  THE  FORM  FACTOR  EQUALS  ZERO 


sin£  =  0. 

sirup  =  0. 

0j  =  S7.9°. 

0j  =  66.4°. 

62  =  S5.S°. 

0,  =  36.9°. 

03  =  83.7°. 

undefined. 

04  =  8t.6°. 

undefined. 

05  =  "9.5°. 

undefined. 

undefined. 

0,  =  75.2°. 

undefined. 
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Figure  3.14  Energy  per  Unit  Solid  Angle,  499-501  MI  lz.,  (Trapezoidal  Form  Factor). 
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IV.  DISCUSSION 


Figures  3.6  through  3.8  show  that  the  energy  radiated  per  unit  solid  angle  from  a 
charge  distribution  with  unity  form  factor  tends  to  be  concentrated  into  an  ever 
narrowing  direction  as  the  range  of  integration  over  frequency  is  increased.  In  the  limit 
of  all  frequencies,  the  radiation  is  confined  to  the  Cerenkov  cone.  The  effect  of 
including  all  frequencies  is  the  same  as  allowing  the  beam  path  length  to  go  to  infinity, 
and  the  resulting  radiation  patterns  are  characteristic  of  that  from  a  point  charge. 

The  zeroes  in  the  radiation  patterns  for  trapezoidal  charge  distributions,  Figures 
3.9  through  3.15,  are  tabulated  in  Tables  1  through  4.  Since  the  tables  and  figures  are 
calculated  from  the  same  theory,  they  are,  of  course,  in  agreement.  Tables  2  and  4 
show  the  zeroes  for  the  "carrier"  oscillation  which  arises  from  the  length  of  the  electron 
beam  path  from  source  to  detector.  Tables  1  and  3  show  the  zeroes  in  the  c  and  <p 
oscillations  resulting  from  the  length  (d  +  b)  and  rise  time  (d-b)  respectively  of  an 
individual  charge  bunch.  Both  Tables  show  zeroes  dictated  by  the  length  of  the  bunch 
charge.  Only  two  zeroes  from  the  pulse  rise  time  occur;  both  at  the  higher  frequency 
of  measurement. 

In  the  cases  studied  here,  the  bunch  length  (d  +  b)  equals  16.5  meters  and  the 
pulse  rise  (d-b)  equals  1.5  meters.  At  the  measuring  frequency  of  100  MHz.  the  ratio 
of  bunch  length  to  wavelength  is  5.5;  large  enough  to  resolve  the  details  of  the  bunch 
size.  Even  at  the  higher  frequency  of  500  MHz.  the  ratio  of  pulse  rise  to  wavelength  is 
2.5  and  consequently  the  rise  and  fall  of  the  pulse  have  little  effect  on  the  radiated 
output.  [Ref.  5] 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


After  studying  the  figures  depicting  the  radiated  Cerenkov  energy  in  this  work, 
one  could  initially  deduce  that  the  use  of  narrowband  detectors  would  be  optimum  for 
detecting  radiation  during  any  electron  accelerator  experiment.  As  shown  and 
previously  discussed,  the  information  produced  is  more  detailed  when  narrowband 
frequency  ranges  are  used  compared  to  wider  band  measurements.  In  particular,  if  one 
frequency  band  is  considered  as  a  function  of  angle,  the  calculated  radiation  energy  has 
much  detail  as  displayed  in  Figures  3.12  and  3.13  The  structure  depends  on  both  the 
path  length  and  bunch  length,  and,  in  principle,  may  yield  information  about  these 
parameters.  One  must  also  consider  the  logistics  of  any  experimental  setup  and  any 
engineering  problems  involved  to  yield  the  maximum  information. 

For  the  broadband  output  with  a  unity  formfactor,  the  energy  output  has  a  peak 
at  the  critical  Cerenkov  angle  as  expected.  This  output  is  typical  of  a  point  charge. 
However,  since  the  nominal  electron  accelerator  bunch  charge  output  has  a  finite 
distribution,  there  is  little  information  to  be  extracted  from  this  particular  output. 

When  the  effect  of  the  form  factor  is  included,  there  exists  a  definite  Cerenkov 
radiation  peak  perpendicular  to  the  path  the  electron  bunch  is  travelling.  This 
phenomenon  may  be  significant  with  respect  to  the  optimum  position  at  which  the 
maximum  energy  produced  by  the  radiation  can  be  detected. 

If  a  narrowband  detection  system  is  realistically  feasible,  it  can  be  used  to 
determine  all  nulls  that  arise  in  the  energy  equation  due  to  the  diffraction  function  and 
the  form  factor.  Although  narrowband  detection  will  produce  the  most  detailed 
output,  its  usefulness  may  be  limited  by  the  availability  of  equipment  designed  to 
discriminate  the  output  within  the  desired  frequency  band  limits. 

In  order  to  obtain  significant  information  about  the  details  of  the  charge 
distribution  within  a  bunch,  it  will  be  necessary  to  perform  experiments  at  frequencies 
such  that  the  ratio  of  rise  time  to  frequency  is  at  least  5.  For  real  beams,  where  the 
rise  times  are  of  the  order  of  a  few  nanoseconds,  this  means  that  narrowband 
measurements  must  be  performed  at  microwave  frequencies. 


To  augment  this  research  the  following  is  recommended: 

1.  Analvze  outputs  with  form  factors  other  than  trapezoidal,  (i.e.  rounded, 
triangular  etc.). 

2.  Determine  the  optimum  (hieher)  frequencv  bands  for  various  frequencv 
intervals  at  which  useful  information  can  be  extracted. 
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APPENDIX 

FORTRAN  PROGRAM 


C  PROGRAM  TO  COMPUTE  THE  ENERGY  RADIATED  PER  UNIT  SOLID  ANGLE 
C  OVER  A  SPECIFIC  FREQUENCY  RANGE 
C 

C  *  ERG  ARRAY  SIZE  MUST  CORRESPOND  TO  THE  NUMBER  OF  * 

C  *  INCREMENTS  OF  THETA  * 

C 

DIMENSION  ERG( 1SOO  ) 

C 

C 

REAL  ERG , COSANG ,COSARG .COSGAR .KKONST >DEE .BEE 
10  REAL  FORFAC.EBEAM, LENGTH, KONST ,C, GEE, TOP, BOT 
20  REAL  ATHETA.RADFUN,  PI,  EREST  ,GAFt1A ,  THETA,  IND.MIDNU 
30  REAL  CONST, MU, CUE, TERG,TTERG,DNU,BNU,DEGINC .DATEND 

INTEGER  IFORM,INCREM,OLOOP .D2LOOP , SELECT , NUMBER 
CHARACTERS  RESPON ,ANSM 
C 
C 

C  ***EBEAM=BEAM  ENERGY  (MEV)  j  FORFAC=FORM  FACTOR  IF(K)) 

C  ***ETA=BEAM  LENGTH/WAVELENGTH  »  BETA=VELOCITY/C( FOR  MEDIUM) 

C  ***GAMMA=RATIO  OF  BEAM  ENERGY  ELECTRON  REST  ENERGY 

c  ***0IFFUN*0I FRACTION  FUNCTION  IIIU))  >  DTHETA=THETA( DEGREES ) 

C  ***THETA=THETA(  RADIANS)  >  LENGTH ^LENGTH  FROM  SOURCE) METERS  ) 

C  ***RADFUN=RADIATION  FUNCTION (  MIDNU=BE GINNING  FREQUENCY 

c  ***CUE=BUNCH  CHARGE  IN  COULOMBS lEREST=  ELECTRON  REST  ENERGY 

C  ***ERG= ARRAY  FOR  STORING  VALUES  OF  RADIATED  ENERGY 

C  »**  FOR  CUMULATIVE  ANGLES  THETA 

C  ***DNU= FREQUENCY  INTERVAL  lTERG*RADIATED  ENERGY  FOR  A 

C  ***  SPECIFIC  ANGLE  THETA 

C 
C 

C  *** ASSIGNMENT  OF  CONSTANTS*** 

C0=2 . 997925E8 
PI=3 . 141S9 
EREST= . 5117 
MU=1 . 2566E-6 
C 
C 

PRINT  * , ' THE  INOEX  OF  REFRACTION  (AIR)  USED  IS  1.000268' 

PRINT  *,'00  YOU  WANT  A  DIFFERENT  VALUE?) Y/N)' 

READ  (5,90 IANSH 
90  FORMAT) Al) 

IF  ( ANSH  .EQ. 'N'  )THEN 
IND=1. 000268 
GO  TO  100 
ELSE 

PRINT  », 'SELECT  THE  INDEX  OF  REFRACTION  (REAL)' 

READ  *>INO 
END  IF 

100  CONTINUE 

C 

C  DETERMINING  THE  SPEED  OF  LIGHT  FOR  THE  MEDIUM 
C  =  CO/IND 
C 

PRINT  *,  'ENTER  BUNCH  CHARGE  (REAL)  IN  COULOFSS' 

READ  *,CUE 
260  CONTINUE 
C 

CCC  COMPUTATION  OF  CONSTANT  (CALLED  Q  IN  THESIS) 

C 


TRI00010 


TRI00020 

TRI00030 

TRI00040 

TRI00050 


CONST *MLWC»( CUE **2 .  I/(8.*(PHh*2.  )) 

CALL  EXCMS  (  ' CLRSCRN'  ) 

C 

C 

PRINT  *,' ENTER  THE  DISTANCE  FROM  THE  SOURCE  I  IN  METERS)' 

PRINT  *,  'THE  RECEIVER  HILL  BE  PLACED' 

READ  »,  LENGTH 

PRINT  •  ,  'SELECT  THE  FORM  FACTOR  DESIRED:' 

PRINT  «,'1=UNITY.  2=TRAPEZOID* 

READ  «,  I FORM 

IF  ( IFORM  .EQ.  11THEN 
F0RFAC=1.0 
GO  TO  290 
ELSE 

PRINT  «, 'ENTER  THE  TOP  AND  BASE  VALUES  (IN  NANOSECONDS)  FOR' 

PRINT  «,'FOR  THE  CURRENT  TRAPEZOID  FUNCTION' 

READ  *,TOP,BOT 
C 

C  THE  NEXT  2  LINES  CONVERT  PULSE  TIME  TO  DISTANCE 
BEE  =  C*TOP/2.E9 
DEE  *  C»BOT/2.E9 

ENO  IF 

290  CONTINUE 

CALL  EXCMS  (  'CLRSCRN'  ) 

C 

C  BEAM  ENERGY  SELECTION:USED  TO  CALCULATE  BETA  A  GAhMA 
C 

PRINT  »,  'SELECT  BEAM  ENERGY:' 

PRINT  * . ' Ia50  MEV  i  2*  SELECT  YOUR  OHN* 

READ  *,IBEAM 

IFIIBEAM  .EG.  1  (THEN 
EBEAM»50. 

GO  TO  300 
ELSE 

PRINT  », 'INPUT  THE  DESIRED  ELECTRON  BEAM  ENERGY  IN  MEV' 

PRINT  •.‘(INCLUDE  DECIMAL  POINT)' 

READ  R.EBEAM 
ENO  IF 

300  CONTINUE 

GAft1A*EBEAM/EREST  TRI00100 

BETA»SQRT(1.-(1./GA»#1A**2.  ))  TRI00I10 

CALL  EXCMS  ( 'CLRSCRN'  ) 

C 

c 

PRINT  •,  ‘SELECT  A  FREQUENCY  BAND  TO  BE  SUM1ED  OVER  AND  THE' 

PRINT  •.'INCREMENT  IN  MHZ:  I  *  A  SUM  OF  FREQUENCIES' 

PRINT  », 'FROM  10  MHZ  TO  100  MHZ  IN  10  MHZ  INCREMENTS' 

PRINT  tt.'Z  «  SELECT  YOUR  OW' 

READ  ».  INCREM 

IF  (INCREM  .EQ.  1  )THEN 
BNU-10.E4 
OLOQP»9 
0NU*10.Et 

CALL  EXCMS  ( 'CLRSCRN'  ) 

GO  TO  400 
ELSE 

CALL  EXCMS  ( 'CLRSCRN'  ) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  CALL  THE  SUBROUTINE  USED  TO  DEFINE  FREQUENCY  BAND  OF  INTEREST  C 
C  AND  THE  INCREMENTAL  STEPS  FOR  SUMflNG  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  FREQ(BNU.DNU.OLOOP) 

EM)  IF 

400  CONTINUE 
C 

C  OUTPUT  GRAPH  AXIS  SELECTION 

PRINT  », 'SELECT  OUTPUT  TYPE' 

PRINT  ». ' 1*£NERGY  VS.  FREQUENCY  FOR  A  SELECTED  ANGLE  THETA' 

PRINT  «,'2«ENERGY( INTEGRATED  OVER  FREQ' 'S  OF  INTEREST)  VS.  THETA' 


PRINT  *,’3=F0RMFACT0R( F<K>)  VS.  FREQ.  FOR  A  SELECTED  ANGLE  THETA1 

READ  *, SELECT 

CALL  EXCMS  <  'CLRSCRN'  ) 

IF  (SELECT  .EG.  1  .OR.  SELECT  .EG.  31  THEN 

GO  TO  470 

ELSE 

GO  TO  450 
END  IF 

450  CONTINUE 

PRINT  *,  'SELECT  A  RANGE  FOR  THETA  AND  THE  INCREMENT  AT  WHICH1 

PRINT  *,  'THE  RADIATEO  ENERGY  WILL  BE  CALCULATED1 

PRINT  *,  '1  =  0  TO  ISO  DEGREES  IN  ONE  DEGREE  INCREMENTS' 

PRINT  *,  'Z  -  SELECT  YOUR  OHN ' 

READ  «,  NUMBER 

IF  (NUMBER  .EG.  I  (THEN 
D2LOOP=181 
DTHETA=0.0 
DEGINC=1 
ELSE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  CALL  THE  SUBROUTINE  USED  TO  DEFINE  THE  INCREMENTAL  C 

C  STEPS  OF  THETA  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CALL  DEGREE  ( DTHETA .DEGINC .D2LOOP  ) 

ENO  IF 

465  CONTINUE 

GO  TO  500 

470  PRINT  #, 'SELECT  AN  ANGLE  (THETA)  OF  INTEREST.1 
READ  *, DTHETA 

MIDNU  s  BNU 
TTERG  =0.0 
0 LOOP =0 LOOP* I 
GO  TO  510 
500  CONTINUE 

DO  600  K=1,D2L00P  TRI00I40 

C 

C  THIS  OUTER  LOOP  HILL  BE  USED  ONLY  FOR  OUTPUT  TYPE  2 
TTERG=0. 

MIDNU=  8NU*( DNU/Z. ) 

510  CONTINUE 

THET  A=0THET  A*PI/180 . 

COSANG  =  COS< THETA) 

KONST  =  2 . »SINI THETA )/( < I ./( IND«6ETA  ) ) -COSANG  ) 

ATHETA  =  (  LENGTH*PI/C )*<  ( l./( IND*BETA )  ) -COSANG  ) 

C 

C  THIS  IMJER  LOOP  IS  USED  FOR  ALL  3  OUTPUT  TYPES  BY  INCREMENTING 
C  FREQUENCY  OVER  THE  CHOSEN  RANGE 

DO  550  J=1,OLOOP 

IFIIFORM  .EG.  1)THEN 
GO  TO  515 

ELSE 

COSARG  *  2.»BEE6PI/C 
COSGAR  *  2.*DEE»PI/C 

KKONST  =  ( C**2.  )/(  2 . «(  PI*»2.  )*(  (  DEE*»2.  )-<  BEE*»2.  ) ) ) 
FORFAC«KKONST»( l./l MIDNU**2 .  )  )*( 1 ./( ABSI COSANG  )»*2 .  )  )* 

D  ( COSI COSARG»CQSANG*MIDNU  l-COSI COSGAR*COSANG*MIDNU  )  ) 

END  IF 

515  CONTINUE 

RADFUN«ABS( KONST»SIN( ATHETARMIDNU )*FORFAC )  TRI001S0 

TERG»CONST*( RADFUN**2.  ) 

IF  (SELECT  .EG.  1)  THEN 

WRITE) 21)520  I  MIDNU, TERG 
520  FORMAT ( 2E 10. 5) 

ELSE  IF  (SELECT  .EG.  3  )THEN 

WRITE  (21,525)  MIONU,FORFAC 
525  FORMAT! E 10 . 5 ,E10 .4 ) 

ELSE 

TTERG  =  TTERG  ♦  TERG 

530  CONTINUE 
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ENO  IF 

MIONU*  MIDNU+ONU 

550  CONTIfAJE  TRI00200 

IF( SELECT  .EG.  1  .OR.  SELECT  .EQ.  5)  THEN 
GO  TO  600 
ELSE 

GO  TO  555 
ENO  IF 

555  CONTINUE 

ERG( K  IsTTERG*  DNU 
GEE =A8S( KONST  1**2. 

HRITEI 21,5601  OTHETA , ERG! K  )  TRI00190 

560  FORMAT  12E10.5  1 

OTHETA  *  OTHETA  ♦  OEGINC 
600  CONTINUE 

C  THE  DATENO  STATEMENT  IS  USED  TO  ENO  THE  DATA  FILE 
C  (NECESSARY  IF  USING  OISSPLA  GRAPHICS) 

DATEN0=999. 99 

WRITE! 21,660  I  OATENO ,OATENO 
660  FORMAT  <2F10.2) 

PRINT  *,  'DO  YOU  WANT  TO  MAKE  ANOTHER  RUNTIY/N)' 

READ  (5,670  1RESP0N 
670  FORMAT! All 

IF! RESPON.EQ. ‘ Y '  I  THEN 
CALL  EXCMS  (  ’CLRSCRN '  1 
GO  TO  100 
ELSE 

GO  TO  700 
ENO  IF 

700  CONTINUE 

STOP  TRI00Z10 

ENO  TRI 00220 

C 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  FR£4<»FRE«,ONU,OIOOP> 

REAL  BFREQ.EFREQ.DIFF.ONU 
INTEGER  OLOOP 

PRINT*  ,' ENTER  BEGIMIN6  FREQUENCY  IN  MHZ!  INCLUDE  DECIMAL  POINT)' 

READ*  .BFREQ 
BFREQ=BFRE«*1.E6 

PRINT*  , 'ENTER  ENDING  FREQUENCY  IN  MHZ I  INCLUDE  DECIMAL  POINT)' 

READ*  .EFREQ 
EFREQ*EFREQ*1.E6 

PRINT*  , 'ENTER  DESIRED  INCREMENT  IN  MHZ! INCLUDE  DECIMAL  POINT)' 

READ*  , DNU 
DNU*0NU*1.E6 
OIFF^EFREQ-BFREQ 
DLOOP*INT( DIFF/DNU  ) 

RETURN 

ENO 

C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  DEGREE ( BDEG , OEGINC  .D2LOOP  ) 

REAL  DOES, EDEG, OEGINC, OIFFF 

INTEGER  D 2 LOOP 

CALL  EXCMS  ('CLRSCRN') 

PRINT  », 'ENTER  STARTING  ANO  ENOING  ANGLES  OF  INTEREST ( REAL  ) ' 

RE AO  *,BOEG,EOEG 

PRINT  »,  'ENTER  THE  INCREMENTAL  VALUE  OF  THETA! DEGREES )  DESIRED.’ 

READ  *,  OEGINC 

OIFFF  »  EDEG-BOEG 

02LOOP  *INT<  DIFFF/OEGINC )»1 

RETURN 

ENO 
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